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MATRIX- VALUED ORTHOGONAL POLYNOMIALS RELATED TO 

(SU(2) x SU(2),diag), II 

ERIK KOELINK, MAARTEN VAN PRUIJSSEN, PABLO ROMAN 

Abstract. In a previous paper we have introduced matrix- valued analogues of the Cheby- 
shcv polynomials by studying matrix- valued spherical functions on SU(2) x SU(2). In par- 
ticular the matrix-size of the polynomials is arbitrarily large. The matrix-valued orthogonal 
polynomials and the corresponding weight function are studied. In particular, we calculate 
the LDU-decomposition of the weight where the matrix entries of L are given in terms of 
Gegenbauer polynomials. The monic matrix- valued orthogonal polynomials P„ are expressed 
in terms of Tirao's matrix- valued hypergeometric function using the matrix- valued differential 
operator of first and second order to which the P„'s are eigenf unctions. From this result we 
obtain an explicit formula for coefficients in the three-term recurrence relation satisfied by the 

w-v | polynomials P n . These differential operators are also crucial in expressing the matrix entries 

of P n L as a product of a Racah and a Gegenbauer polynomial. We also present a group 
theoretic derivation of the matrix-valued differential operators by considering the Casimir 

^ ■ operators corresponding to SU(2) x SU(2). 

U: 

!: 

_., ! 1. Introduction 

Matrix-valued orthogonal polynomials have been studied from different perspectives in re- 
■^- ■ cent years. Originally they have been introduced by Krein [T7j, [18]. Matrix- valued orthogonal 

polynomials have been related to various different subjects, such as higher-order recurrence 
equations, spectral decompositions, and representation theory. The matrix-valued orthog- 
onal polynomials studied in this paper arise from the representation theory of the group 
SU(2) x SU(2) with the compact subgroup SU(2) embedded diagonally, see [15] for this par- 
ticular case and Gangolli and Varadarajan [5], Tirao [20], Warner [22J for general group 
theoretic interpretations of matrix-valued spherical functions. An important example is the 
study of the matrix- valued orthogonal polynomials for the case (SU(3), U(2)), which has been 
studied by Griinbaum, Pacharoni and Tirao [H] mainly exploiting the invariant differential 
operators. In [15J we have studied the matrix-valued orthogonal operators related to the case 
(SU(2) x SU(2), SU(2)), which lead to the matrix-valued orthogonal polynomial analogues of 
Chebyshev polynomials of the second kind U n , in a different fashion. In the current paper we 
study these matrix-valued orthogonal polynomials in more detail. 
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In order to state the most important results for these matrix- valued orthogonal polynomials 
we recall the weight function [T5J Thm. 5.4]: 

m 

W{x) n . m = Vl - x 2 y^ a t (m, n) U n+m _ 2t (x), 

i=0 (1.1) 

. {2£+l){2£-m)\m\, ^ m _ t (n - 2£) m . t (2£ + 2 - t) t 

atM = —l W)^ { ~ 1) (n + 2) m . t tl 

if n > m and W(x) nym = W{x) m , n otherwise. Here and elsewhere in this paper £ G |N, 
n, m G {0, 1, • • • , 2£}, and U n is the Chebyshev polynomial of the second kind. Note that the 
sum in (II. ip actually starts at min(0,n + m — 2£). It follows that W : [— 1, 1] — > M2e+i(C), 
W(x) = \W(x) n ,m) nm=0 i is a (2£ + l) x (2£ + l)-matrix- valued integrable function such that all 

moments j_ 1 x n W (x) dx, n£N, exist. From the construction given in [T5j §5] it follows W(x) 
is positive definite almost everywhere. By general considerations, e.g. [TU], we can construct 
the corresponding monic matrix-valued orthogonal polynomials {P n }^ =0 , so 

(P n , P m ) w = f P n {x) W(x) (P m {x))*dx = S nm H n , 0<H n e M 2m (C) (1.2) 

where H n > means that H n is a positive definite matrix, P n (x) = Yl^o^-^'k with P£ e 
M 2 ^_)_i(C) and P™ = I, the identity matrix. The polynomials P n are the monic variants of 
the matrix-valued orthogonal polynomials constructed in [T5] from representation theoretic 
considerations. Note that ( II. 2p defines a matrix-valued inner product (•, -)w on the matrix- 
valued polynomials. Using the orthogonality relations for the Chebyshev polynomials U n it 
follows that 

(H ) nm = ^2(n + l)(2£-n + l) (L3) 

which is in accordance with [T5J Prop. 4.6]. From [T5] we can also obtain an expression for 
H n by translating the result of [T5, Prop. 4.6] to the monic case in [T5J (4.6)], but since the 
matrix T^ in [TBI (4.6)] is relatively complicated this leads to a complicated expression for the 
squared norm matrix H n in (II. 2p . In Corollary 15.41 we give a simpler expression for H n from 
the three-term recurrence relation. 

These polynomials have a group theoretic interpretation as matrix-valued spherical func- 
tions associated to (SU(2) x SU(2), SU(2)), see [15] and Section [71 In particular, in [151 §5] 
we have shown that the corresponding orthogonal polynomials are not irreducible, but can 
be written as a 2-block-diagonal matrix of irreducible matrix- valued orthogonal polynomials. 
Indeed, if we put J G M^(C), J nm = <5 n +m,2£ we have JW(x) = W(x)J for all x G [— 1, 1], and 
by [13 Prop. 5.5] J and I span the commutant {Y G M 2£ (C) | [Y, W{x)] = OVa; G [-1,1]}. 
Note that J is a self-adjoint involution, J 2 = I, J* = J. It is easier to study the polynomials 
P n , and we discuss the relation to the irreducible cases when appropriate. 

In this paper we continue the study of the matrix-valued orthogonal polynomials and the 
related weight function. Let us discuss in some more detail the results we obtain in this paper. 
Some of these results are obtained employing the group theoretic interpretation and some are 
obtained using special functions. Essentially, we obtain the following results for the weight 
function: 
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(a) explicit expression for det(W(x)), hence proving [151 Conjecture 5.8], see Corollary 12. 3 j 

(b) an LDU-decomposition for W in terms of Gegenbauer polynomials, see Theorem 12.11 

Part (a) can be proved by a group theoretic consideration, and gives an alternative proof 
for a related statement by Koornwinder [16], but we actually calculate it directly from (b). 
The LDU-decomposition hinges on expressing the integral of the product of two Gegenbauer 
polynomials and a Chebyshev polynomial as a Racah polynomial, see Lemma 12.71 
For the matrix-valued orthogonal polynomials we obtain the following results: 

(i) P n as eigenfunctions to a second-order matrix-valued differential operator D and a first- 
order matrix- valued differential operator E, compare [T5j §7], see Theorem 13 . 1 1 and Sec- 
tion H2 

(ii) the group-theoretic interpretation of D and E using the Casimir operators for SU(2) x 
SU(2), see Section [TJ for which the paper by Casselman and Milicic [5] is essential; 

(iii) explicit expressions for the matrix entries of the polynomials P n in terms of matrix- valued 
hypergeometric series using the matrix- valued differential operators, see Theorem 14. 5} 

(iv) explicit expressions for the matrix entries of the polynomials P n L in terms of (scalar- 
valued) Gegenbauer polynomials and Racah polynomials using the LDU-decomposition 
of the weight W and differential operators, see Theorem 16. 2( 
(v) explicit expression for the three-term recurrence satisfied by P n , see Theorem 15.31 

In particular, (i) and (ii) follow from group theoretic considerations, see Section |3] and 
This then gives the opportunity to link the polynomials to the matrix-valued hypergeometric 
differential operator, leading to (iii). The explicit expression in (iv) involving Gegenbauer 
polynomials is obtained by using the LDU-decomposition of the weight matrix and the differ- 
ential operator D. The expression of the coefficients as Racah polynomials involves the first 
order differential operator as well. Finally, in [15, Thm. 4.8] we have obtained an expression 
for the coefficients of the three-term recurrence relation where the matrix entries of the coeffi- 
cient matrices are given as sums of products of Clebsch-Gordan coefficients, and the purpose 
of (v) is to give a closed expression for these matrices. The case £ = 0, or the spherical 
case, corresponds to the Chebyshev polynomials U n (x), which occur as spherical functions 
for (SU(2) x SU(2), SU(2)) or equivalently as characters on SU(2). For these cases almost all 
of the statements above reduce to well-known statements for Chebyshev polynomials, except 
that the first order differential has no meaning for this special case. 

The structure of the paper is as follows. In Section [2] we discuss the LDU-decomposition 
of the weight, but the main core of the proof is referred to |A] In Section [3] we discuss the 
matrix-valued differential operators to which the matrix-valued orthogonal polynomials are 
eigenfunctions. We give a group theoretic proof of this result in Section [71 In [T51 §7] we 
have derived the same operators by a judicious guess and next proving the result. In order to 
connect to Tirao's matrix-valued hypergeometric series, we switch to another variable. The 
connection is made precise in Section HJ This result is next used in Section [5] to derive a 
simple expression for the coefficients in the three-term recurrence of the monic orthogonal 
polynomials, improving a lot on the corresponding result [151 Thm. 4.8]. In Section [6] we 
explicitly establish that the entries of the matrix-valued orthogonal polynomials times the 
L-part of the LDU-decomposition of the weight W can be given explicitly as a product of 
a Racah polynomial and a Gegenbauer polynomial, see Theorem 16.21 Some of the above 
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statements require somewhat lengthy and/or tedious manipulations, and in order to deal with 
these computations and also for various other checks we have used computer algebra. 

As mentioned before, we consider the matrix-valued orthogonal polynomials studied in this 
paper as matrix- valued analogues of the Chebyshev polynomials of the second kind. As is well 
known, the group theoretic interpretation of the Chebyshev polynomials, or more generally 
of spherical functions, leads to more information on these special functions, and it remains to 
study which of these properties can be extended in this way to the explicit set of matrix- valued 
orthogonal polynomials studied in this paper. This paper is mainly analytic in nature, and 
we only use the group theoretic interpretation to give a new way on how to obtain the first 
and second order matrix- valued differential operator which have the matrix- valued orthogonal 
polynomials as eigenfunctions. We note that all differential operators act on the right. The 
fact that we have both a first and a second order differential operator makes it possible to 
consider linear combinations, and this is useful in Section H] to link to Tirao's matrix-valued 
differential hypergeometric function and Section [6] in order to diagonalise (or decouple) a 
suitable matrix- valued differential differential operator. 

We finally remark that J. A. Tirao has informed us that Ignacio Zurrian has obtained results 
of a similar nature by considering matrix- valued orthogonal polynomials for the closely related 
pair (SO(4), SO(3)). We stress that our results and the results by Zurrian have been obtained 
independently. 



2. LDU-DECOMPOSITION OF THE WEIGHT 

In this section we state the LDU-decomposition of the weight matrix W in (11.11) is discussed. 
The details of the proof, involving summation and transformation formulas for hypergeometric 
series (up to 7_F 6 -level), is presented in|A] Some direct consequences of the LDU-decomposition 
are discussed. The explicit decomposition is a crucial ingredient in Section[6l where the matrix- 
valued orthogonal polynomials are related to the classical Gegenbauer and Racah polynomials. 

In order to formulate the result we need the Gegenbauer, or ultraspherical, polynomials, 
see e.g. [2], [II], [H], defined by 

Cr'M^^f-"-'^ 2 ";!^). (2.1) 

The Gegenbauer polynomials are orthogonal polynomials; 



(1 x ) C n (x)C m (x)dx - o nm nl{n + a)T{a) 

nT(n + 2a)2 1 - 2a 
r(a) 2 {n + a) n\ 

Theorem 2.1. The weight matrix W has the following LDU-decomposition; 



(2.2) 



W(x) = Vl - x 2 L(x) T{x) L(xY, are [-1,1], 
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where L: [— 1, 1] — > M2£+i(C) is the unipotent lower triangular matrix 

0, k > m 



L(x) mk ={ m!(2fc + l)! 



C^lHx), k<m 



and T: [—1,1] — > M2£+i(C) is the diagonal matrix 

T(x) kk = cM(l - x'f cM = 4»(H)'(2* + D(2< + *+l)l(2<-*)l 
1 '" t{,( >' k( > ((2* + 1)1)2 ((2«)!)2 

Note that the matrix-entries of L are independent of £, hence of the size of the matrix- valued 
weight W. Using 

d k 



dx k 



C^{x) = 2\ a ) k C, t 



k („,\ n( a+k \ 



n—k 



X) 



we can write uniformly L(x) m k = rrks (m+k+iv. ^d^" ! 3 -)- In Theorem 16.21 we extend Theorem 
I2.1[ but Theorem 12.11 is an essential ingredient in Theorem 16.21 

Since W(x) is symmetric, it suffices to consider (n, m)-matrix-entry for m < n of Theorem 
12.11 Hence Theorem 12.11 follows directly from Proposition 12.21 using the explicit expression 
(II. ip for the weight W. 

Proposition 2.2. The following relation 

m m 

£ a t (m, n)U n+m . 2t {x) = £ (3 k (m, n)(l - x 2 ) fc C^^C^ix) 

t=0 fc=0 

with the coefficients at(m,n) given by (II. ip and 

o< ml nl ,,,.o2fc/o, ,A2£ + k + l)\ (2£-k)\ 

8 k (m,n) = -, : -r- : —A; A; 2 2fe (2A; + 1)- , ', \, - 

' K ' ; (m + k + l)\(n + k + l)\ v ; {2£)\{2£)\ 

holds for all integers < m < n < 2£, and all £ G |N. 



Before discussing the proof we list some corollaries of Theorem 12.11 First of all, we can use 
Theorem 12. II to prove [15], Conjecture 5.8], see (a) of Section [TJ 

Corollary 2.3. det(W{x)) = (1 - x 2 ) 2 ^) 2 Y^ =Q c k {£). 

Remark 2.4. We also have another proof of this fact using a group theoretic approach to 
calculate det($o(aO), see [15] and Section [7] for the definition of $o, and W is up to trivial 
factors equal to ($o)( < l ) o)*- This proof is along the lines of Koornwinder |16| . 

Secondly, using J e M 2 £+i(C), J nm = 5 n +m,2e. and W(x) = JW(x)J, see [El Prop. 5.5, 
§6.2], we obtain from Theorem 12. II the UDL-decomposition for W. For later reference we also 
recall JP n (x)J = P n (x), since both are the monic matrix- valued orthogonal polynomials with 
respect to W(x) = JW(x)J. 



Corollary 2.5. W{x) = y/l - x 2 (JL(x)J) (JT(x)J) (JL(x)J) , x e [-1,1] gives the UDL- 
decomposition of the weight W . 
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Thirdly, considering the Fourier expansion of the weight function W(cos6), and using the 
expression of the weight in terms of Clebsch-Gordan coefficients, see [151 (5-4), (5.6), (5.7)] 
we obtain a Fourier expansion, which is actually equivalent to Theorem 12.11 

Corollary 2.6. We have the following Fourier expansion 

yv i)k(2k , (m-k + l) k (n-k + l) k (2e+k+l)\(2e-k)\ t(n+m)t 
£-T (m+l) k+1 (n + l) k+1 (2£)\(2£)\ 



2it\2k 77 I* n '^ + l.„2Jt\ 77 I * m 5^ + l.„2it 



fc=0 

\ — n / \ — m 

21 n 2l-n m 2£-m (n\/2t-n\ tm\ f2i-m\ 
V^V^ V^ V^ V Wj ^ 2 1 Mil!: ^ i ^((n-ji+J2)-(m-ii+z 2 ))t 

A^ A^ A^ A-^ Z^ pA pA e 

j=0ji=0 J2=0 U=0 12=0 \j/ VjrV 

jl+J2=j il+l2=j 

Proof. In [T5l §5, 6] the weight function Vy(cost) was initially defined as a Fourier polynomial 
with the coefficients given in terms of Clebsch-Gordan coefficients. After relabeling this gives 

21 n 2£-n m 2l-m (n\ (2l-n\ (m\ (2i-m 



EEEEE ■ 



tn\(2£-n\ /m\ f2H.-m\ 

W V 31 I W V '2 ' r i((n-ji+i 2 )-(m-i 1 +i 2 ))t 



j=0 jl=0 J2=0 U=0 «2=0 \j / \j 

3\+32=3 U+J2=j 

min(m,n) 

(L(cost)T(cost)L(cost)') nm = £ &(ro,n) sin 2fc t cf^ (cos t)tf£+i> (cost) 

fc=0 

where we have used [T5J (5.10)] to express the Clebsch-Gordan coefficients in terms of binomial 
coefficients. 

Using the result [31 Cor. 6.3] by Koornwinder and Badertscher together with the Fourier 
expansion of the Gegenbauer polynomial, see [H (2.8)], [21 (6.4.11)], [HI (4.5.13)], we find the 
Fourier expansion of sin fc t C^_ k (cos t) in terms of Hahn polynomials defined by 

Q k U;a,p,N)= 3 F 2 (~ Kk+ ^^ h ~ J ;l), ke{0,l,---,N}, (2.3) 

V a + l,-N J 

see [21 p. 345], [EEJ §6.2], pH §1.5]. For A = 1 the explicit formula is 



i k {n + l) fc+1 {n-k)\ k (fc+1) _ ^ nk*^'-^* 

— n 



(2.4) 



using the generating function [14, (1.6.12)] for the Hahn polynomials in the last equality. 
Plugging this in the identity gives the required result. □ 

In the proof of Proposition 12.21 and Theorem 12.11 given in |A] we use a somewhat unusual 
integral representation of a Racah polynomial. Recall the Racah polynomials, [21 p. 344], [131 
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§1.2], defined by 

p f\M R x\ f f-k,k + a + B + l,-t,t + 7 + S + l \ 

R k A (t );a,8,j,S = 4^3 .ifl.x.i , 1 ;1 2 " 5 

\ a + l,p + + 1,7 + 1 / 

where \(t) = t(t + 7 + S + 1), and one out of a + 1, /3 + <5+l, 7 + I equals — iV with a 
non-negative integer iV. The Racah polynomials with < k < N form a set of orthogonal 
polynomials for t G {0, 1, ■ ■ • , iV} for suitable conditions on the parameters. For the special 
case of the Racah polynomials in Lemma 12.71 the orthogonality relations are given in 1X1 

Lemma 2.7. For integers 0<t,k<m<nwe have 

1 n _ r 2 ) k+ '2r (k+1) (r)r {k+1) (r)n m Ht - v /?r ( /c + D (k + i) m - k 

(A; + l)„_ fc (-l) fc (2A; + 2) m+n _ 2fc (A; + 1)! 

x ^ — ir — \ , 7?rr -Rk(Kt); o, o, -n - 1, -m - 1 

(n — k)\ (n + m + 1)! 



Remark 2.8. Lemma 12. 71 can be extended using the same method of proof to 

1 (i - x 2 r +k ^ctT 1 \ x )ctT 1 \ x )ci p l m ^ t {x) dx = 



-1 



(a + k + l) m - fc (2fc + 2a + 2) n _ fc (-m + g - a - l) m _ t (8) n ^ t ^T(a + k + |) 
(m — A;)! (n — m)! {m — t)\ T{ot + n + m — t + 2) 

/A; — ?ti, — m — 2a — k — l,t — m, 3 + n — t 
\ p — a — 1 — m, —m — a, n — m + 1 



(2.6) 



assuming n > m. Lemma [2.71 corresponds to the case a = 0, 8 = 1 after using a transforma- 
tion for a balanced 4i*3-series. Note that 4i*3-series can be expressed as a Racah polynomial 
orthogonal on {0, 1, • • ■ , m} in case a — 0ov8 — a+1, which corresponds to Lemma 12.71 
We do not use ( 12 .6p in the paper, and a proof follows the lines of the proof of Lemma 12.71 as 
given in LAI 

In [151 Thm. 6.5], see SectionLTl we have proved that the weight function W is not irreducible, 
meaning that there exists Y G M 2 ^+i(C) so that 

YW{x)Y* = (^ w ° 2{x ^j , YY* = I = Y*Y (2.7) 

and that there is no further reduction. 
Writing 

Y ~ ( A B \ L(x) - ( h{x) ° ^1 T(x) - ( h{x) ° 
" \C D) ' MXj ~ \r(x) h{x)) ' K ' ' \ t 2 {x) 

with A, D diagonal and B, C antidiagonal, see [T5j Cor. 5.6], h{x), l 2 (x) lower-diagonal matri- 
ces and r(x) a full matrix, we can work out the block- diagonal structure of YL(x)T(x)L(x) t Y t . 
It follows that the off-diagonal blocks being zero is equivalent to 

(AZi(a;)ti(aO + Sr(x)ti(x)) {l^xfC 1 + r(aO*Z>*) + Bl 2 (x)t 2 (x)l 2 (xY D l = 0. (2.8) 
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This can be rewritten as an identity for four sums of products of two Gegenbauer polynomials 
involving the weight function and the constants in Theorem 12. II being zero. We do not write 
the explicit results, since we do not need them. 



3. Matrix- valued orthogonal polynomials as eigenfunctions of 
matrix- valued differential operators 

In [151 §7] we have derived that the matrix-valued orthogonal polynomials are eigenfunc- 
tions for a second and a first order matrix-valued differential operator by looking for suitable 
matrix-valued differential operators self-adjoint with respect to the matrix-valued inner prod- 
uct (-, -)\y. The method was to establish relations between the coefficients of the differential 
operators and the weight W, next judiciously guessing the general result and next proving it 
by a verification. In this paper we show that essentially these operators can be obtained from 
the group theoretic interpretation by establishing that the matrix-valued differential operators 
are obtainable from the Casimir operators for SU(2) x SU(2). Since the paper is split into a 
first part of analytic nature and a second part of group theoretic nature, we state the result 
in this section whereas the proofs are given in Section [JJ Sections H] and |6] depend strongly on 
the matrix- valued differential operators in Theorem 13.11 

Recall that all differential operators act on the right, so for a matrix- valued polynomial 
P: M — y Mjy(C) depending on the variable x, the s-th order differential operator D = 
E!=o 7&Fi(?), Fi-.R^ M N (C), acts by 

(PD)(x) = J2^(x)F i (x), PD:R^M N (C) 

i=0 

where (^r(%)) — rf '^ m (x) is a matrix which is multiplied from the right by the matrix Fi(x). 
The matrix-valued orthogonal polynomial is an eigenfunction of a matrix-valued differential 
operator if there exists a matrix A 6 Mjv(C), the eigenvalue matrix, so that PD = AP as 
matrix- valued functions. Note that the eigenvalue matrix is multiplied from the left. For more 
information on differential operators for matrix- valued functions, see e.g. [10], [21] . 

We denote by E^ the standard matrix units, i.e. Eij is the matrix with all matrix entries 
equal to zero, except for the (i,j)-ih entry which is 1. By convention, if either i or j is not in 
the appropriate range, the matrix E^ is zero. 

Theorem 3.1. Define the second order matrix-valued differential operator 

2> d 2 ( d 



2£ 21 21 

C = Y,(M- i)Ei, i+ i + J2 iE ^ U = (2£ + 3)1, V = -J2 1(2£ - i)E v . 

i=0 i=0 i=0 
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and the first order matrix-valued differential operators 

*=(S ( * o+xJi)+i 

~_ ^ (2£ -i) ^ (£-z) V i 

i=0 " i=0 i=0 

^i - - 2^ ~~f~ tji ^ A ~l^ z^ hi » 

j=0 i=0 

t/iera i/ie monic orthogonal matrix-valued orthogonal polynomials P n satisfy 

2£ 

P n D = K{D)P n , An{D) = Y, (" n ( n " X ) - n ( 2i + 3 ) + *( 2£ - *)) E ih 

i=0 

r, * a ,^„ a /^ ^fnU-i) (2£ + 2)(i - 2£)\ „ 
P n £ = A n (P)P n , A n (P) = ^ ( -A— -2 - * -A i #... 

«=o ^ ' 



and £/te operators D and E commute. The operators are symmetric with respect to W . 

The group theoretic proof of Theorem 13.11 is given in Section [71 Theorem 13.11 has been 
proved in [151 Thms. 7.5, 7.6] analytically. The symmetry of the operators with respect to 
W means that (PD,Q)w = {P,QD)w and {PE,Q)w = {P,QE)w for all matrix-valued 
polynomials with respect to the matrix- valued inner product (-,-)w defined in (\1.2\i . The last 
statement follows immediately from the first by the results of Griinbaum and Tirao JTU] . Also, 
[D, E] = follows from the fact that the eigenvalue matrices commute. In the notation of [TU] 
we have D,E e T>(W), where T>{W) is the *-algebra of matrix- valued differential operators 
having the matrix-valued orthogonal polynomials as eigenfunctions. 

Note that E has no analogue in case £ = 0, whereas D reduces to the hypergeometric 
differential operator for the Chebyshev polynomials U n . 

The matrix-differential operator D is J-invariant, i.e. JDJ = D. The operator E is almost 
J- anti- invariant, up to a multiple of the identity. This is explained in Theorem 17.151 and the 
discussion following this theorem. In particular, D descends to the corresponding irreducible 
matrix- valued orthogonal polynomials, but E does not, see also [151 §7]. 

4. Matrix- valued orthogonal polynomials as matrix- valued hypergeometric 

functions 

The polynomial solutions to the hypergeometric differential equation, see f)6.5p . are uniquely 
determined. Many classical orthogonal polynomials, such as the Jacobi, Hermite, Laguerre 
and Chebyshev, can be written in terms of hypergeometric series. For matrix-valued valued 
functions Tirao [21] has introduced a matrix-valued hypergeometric differential operator and 
its solutions. The purpose of this section is to link the monic matrix-valued orthogonal 
polynomials to Tirao's matrix-valued hypergeometric functions. 
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We want to use Theorem 13.11 in order to express the matrix- valued orthogonal polynomials 
as matrix- valued hypergeometric functions using Tirao's approach [21]. In order to do so we 
have to switch from the interval [— 1, 1] to [0, 1] using x = 1 — 2u. We define 

R n (u) = (-1)"2-"P„(1 - 2u), Z(u) = W{1 - 2«) (4.1) 

so that the rescaled monic matrix-valued orthogonal polynomials R n satisfy 

Z(u) = W(l-2u), [ R n {v)Z{u)R m {u)*du = 2- l - 2n H n . (4.2) 



In the remainder of Section |4] we work with the polynomials R n on the interval [0, 1]. It is a 
straightforward check to rewrite Theorem 13.11 

Corollary 4.1. Let D and E be the matrix valued differential operators 

D = ^~ «>;£ + (I) < c - Mt/ » " v ' E = (1) <" B ' + B °> + A - 

where the matrices C, U , V , B , B\ and A are given by 



^ (2t-i) ^,(2£ + 3) * 

2^ 2 ' m ^ 2 ~ ^ 

i=0 i=0 i=0 

21 21 , n „ . _ w . „«> 2^ 



j=0 i=0 i=0 

i=0 ' i=0 i=0 

JTien D and E are symmetric with respect to the weight W , and D and E commute. Moreover 
for every integer n > 0, 

R„D = A n {D)R n , A n {D) = £ (-n(n - 1) - n(2£ + 3) + i(2t - %)) E it , 

i=0 

n r, » ,^„ » ,^ v^/ nU-i) (2£ + 2)(i - 2£)\ „ 
i? n £ = A n {E)R n , A n (E) = J2 I V^ + " ^ " J #«■ 

i=0 ^ ' 

It turns out that it is more convenient to work with D a = D + aE for a 6 R, so that 
R n D a = A n (D a )R n with diagonal eigenvalue matrix A n (D a ) = A n (D) + aA n (E). By [10, 
Prop. 2.6] we have A n (D a ) = —n 2 — n(U a — 1) — V a . Since the eigenvalue matrix A n (D a ) is 
diagonal, the matrix-valued differential equation R n D a = A n (D a )R n can be read as 21 + 1 
differential equations for the rows of R n . The z-th row of R n is a solution to 

u(l-u)p"(u)+p'(u)(C a -uU a )-p(u)(V a + X) = 0, X= (A n (D a )).. (4.3) 

for p: C — > C 2£+1 a (row-)vector- valued polynomial function. Here C a = C + tt-Bo? K* = 
£/ — aBi, V a = V — aA using the notation of Corollary 14. 1[ Now (14. 3 j) allows us to connect 
to Tirao's matrix- valued hypergeometric function [2T], which we briefly recall in Remark 14.21 
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Remark 4.2. Given d x d matrices C, U and V we can consider the differential equation 

z(l - z)F"{z) + (C - zU)F'(z) - VF(z) = 0, z e C, (4.4) 

where F: C — > C d is a (column-)vector- valued function which is twice different iable. It is 
shown by Tirao [2T] that if the eigenvalues of C are not in — N, then the matrix valued 
hypergeometric function 2 H\ defined as the power series 

oo 



\ / i— n 



(4-5) 

[C,U,V] = 1, [C,U,V] i+1 = (C + i)- 1 {i 2 + i(U -l) + V)[C,U,V\i 
converges for \z\ < 1 in M<j(C). Moreover, for Fq G C d the (column-)vector- valued function 

is a solution to (14.4p which is analytic for \z\ < 1, and any analytic (on \z\ < 1) solution to 
(14. 4 p is of this form. 

Comparing Tirao's matrix-valued hypergeometric differential equation (14.41) with (14.31) and 
using Remark I4.2[ we see that 

Po £ C 2 ^ +1 , are the solutions to (14. 3 p which are analytic in |w| < 1 assuming that eigenvalues 
of C* are not in — N. We first verify this assumption. Even though V a and U a are symmetric, 
we keep the notation for transposed matrices for notational esthetics. 

Lemma 4.3. For every £ G |N, the matrix C a is a diagonalizable matrix with eigenvalues 
(2j + 3)/2,je{0,...,2£}. 

Proof. Note that C a is tridiagonal, so that V\ = J2 n =oPnW e n is an eigenvector for C a for the 
eigenvalue A if and only if 

3 . (2£-n)(2£ + a) 

-(A--)Pn(A) = ^ Pn+l(A) 

(2£ + a)(2£-n) + n(2£-a)\ ,,, (2£ - a)n 

~ ^ " "J P-W + U P^(\). 

The three-term recurrence relation corresponds precisely to the three-term recurrence relation 
for the Krawtchouk polynomials for iVeN, 

_' N ;-), n,xG{o,i,---,iV}, 

see e.g. p. 347], HU §6.2], [H §1.10], with N = 2£, p = 2£ts. The Krawtchouk polynomials 

are orthogonal with respect to the binomial distribution for < p< 1, or a G (— 2£, 2£), and 

we find 

2/4- rv 
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and the eigenvalues of C a are | 4- x, x G {0, 1, • • • ,2£}. This proves the statement for 
ae(-2£,2£). 

Note that for a ^ ±2£, the matrix C a is tridiagonal, and the eigenvalue equation is solved 
by the same contiguous relation for the 2-^1-series leading to the same statement for \a\ > 21. 
In case a = ±2£ the matrix C±2i is upper or lower triangular, and the eigenvalues can be read 
off from the diagonal. □ 

In particular, we can give the eigenvectors of C a explicitly in terms of terminating 2-F1- 
hypergeometric series, but we do not use the result in the paper. 

So ( 14. 6 p is valid and this gives a series representation for the rows of the monic polynomial 
R n . Since each row is polynomial, the series has to terminate. This implies that there exists 
n G N so that [C*, [/*, V l a + A] n+1 is singular and ^ P G Ker ([C*, [/*, V l a + A] n+1 ). 

Suppose that n is the least integer for which [C^, £/*, V^ + A] n+ i is singular, i.e. [C*, [/*, V* + 
\}i is regular for all i < n. Since 

[Cl Ul Vi + A] n+1 = (C* + n)- 1 (n 2 + n(Ui - 1) + V* + A) [C£, E& ^ + A] n (4.7) 

and since the matrix {C a + n) is invertible by Lemma I4.3[ [C l a) U^,V£ + A] n+ i is a singular 
matrix if and only if the diagonal matrix 

M°(A) = (n 2 + n(Ui - 1) + ^ + A) 

= (n 2 + n([/ Q -l) + \/ a + A) = A-A n (D Q ) 

is singular. Note that the diagonal entries of M°(A) are of the form A — A°(n), so that M n (X) 
is singular if and only if A = \f(n) for some j G {0, 1, ■ • • , 2£}. We need that the eigenvalues 
are sufficiently generic. 

Lemma 4.4. Let a G R \ Q. Then (j, n) — (i, to) G {0, 1 • • • , 2£} x N z'/ and only if AJ(n) = 
A?(m). 



- A?(n) - Af (to) - (m - n)(n + m + 2 + 2£) + (j - z) I °^ * ^ - j - i + 2£ 



Proof. Assume A"(n) = Xfijn) and let (j, n), (z, m) G {0, 1 • • • , 2£} x N, then 

A?(Wl = fro, - n\(n 4- m 4- 5> 4- 9/"l 4- fi - iW — 

2£ 

If j 7^ z, then we solve for a = (2 e+4) ( ~ _^f + + + j + i — 2£\ which is rational. 

Assume next that j = i, then (m — n) (n 4- m + 2 + 2£) = 0. Since n,m,£ > 0, it follows 
that n = m and hence (j, n) = (i,m). □ 

Assume a irrational, so that Lemma l4~4l shows that M n (A"(ra)) is singular if and only if n = 
m. So in the series f 14 . 6 1) the matrix [C£, £/* , V^ + A] n +i is singular and [C* , [/*, V^ + A]j is non- 
singular for < z < n. Furthermore, by Lemma 14741 we see that the kernel of [C£, L/£, V^+A] n+ i 
is one-dimensional if and only if A = A"(z), z G {0, 1, • ■ ■ , 2£}. In case A = A"(z) 

P = [ClU t a X + KWn 1 e l 

is determined uniquely up to a scalar, where e* is the standard basis vector. 

We can now state the main result of this section, expressing the monic polynomials R n as 
a matrix-valued hypergeometric function. 
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Theorem 4.5. With the notation of Remark \4.2\ the monic matrix-valued orthogonal polyno- 
mials can be written as 

(Rn{u)) i3 = (aHi^'^^^J^nia^-Kl + W)]^)' 

for all a6i 

Note that the left hand side is independent of a, which is not obvious for the right hand 
side. 

Proof. Let us first assume that a is irrational, so that the result follows from the considerations 
in this section using that the z-th row of R n {u) is a polynomial of (precise) degree n. The 
constant follows from monocity of R n , so that (R n (u))a = u n . 

Note that the left hand side is independent of a, and the right hand side is continuous in 
a. Hence the result follows for afl. □ 

In the scalar case £ = Theorem 14.51 reduces to 

R n {u) = (-A)- n (n+l) 2 F 1 (~ n ^ + 2 ;u\ (4.9) 

which is the well-known hypergeometric expression for the monic Chebyshev polynomials, see 
U §2.5], HH (4.5.21)], [14, (1.8.31)]. 

5. Three- term recurrence relation 

Matrix- valued orthogonal polynomials satisfy a three-term recurrence relation, see e.g. [6], 
[10] . In [151 Thm. 4.8] we have determined the three-term recurrence relation for the closely 
related matrix-valued orthogonal polynomials explicitly in terms of Clebsch-Gordan coeffi- 
cients. The matrix entries of the matrices occurring in the three-term recurrence relation 
have been given explicitly as sums of products of Clebsch-Gordan coefficients. The purpose of 
this section is to give simpler expressions for the monic matrix-valued orthogonal polynomials 
using the explicit expression in terms of Tirao's matrix-valued hypergeometric functions as 
established in Theorem 14.51 

From general theory the monic orthogonal polynomials R n : R — > Mjy(C) satisfy a three- 
term recurrence relation uR n {u) = R n+ \{u) + X n R n (u) + Y n R n -i(u), n > 0, where i?_i = 
and X n , Y n e M 2 ^+i(C) are matrices depending on n and not on x. Lemma [54] should be 
compared to P Lemma 2.6]. 

Lemma 5.1. Let {-R n }n>o be the sequence of monic orthogonal polynomials and write R n {u) = 
^fc =0 -Rfc uk ; R-k e M N (C), and R™ = I. Then the coefficients X n , Z n of the three-term 
recurrence relation are given by 

Proof. Let (-, •) denote the matrix- valued inner product for which the monic polynomials are 
orthogonal. Using the three-term recursion, orthogonality relations and expanding the monic 
polynomial of degree n + 1 gives 

(uR n — X n R n — Y n R n _i,R n ) = (R n+ i,R n ) = (u n+ ,R n ) + _R™ + (u n } R n ). 
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By the orthogonality relations the left hand side can be evaluated as 

(uR n — X n R n — Y n Rn-i, R n ) = (u n , R n ) + R^ l _ 1 (u n ,R n ) — X n (R n ,R n ) 

and comparing the two right hand sides gives the required expression for X n , since (R n , R n ) = 
(u n , R n ) is invertible. 
The expression for Y n follows by considering on the one hand 

(uR n — X n R n — Y n R n _i, R n -\) = (Rn+i, R n -i) 

= (u n+1 ,R n . 1 ) + K +1 ( M ",i? n _!) + Rlt\{u n -\R n ^). 

while on the other hand the left hand side also equals 

(u n+ \R n ^) + Rl^^R^) + R n n _ 2 {u n -\R n ^) 

— X n (u n ,R n -i) — X n R™_ l (u n ~ ,R n -i) — Y n (R n _i, R n _i) 

and using the expression for X n and cancelling common terms gives the required expression, 
since (R n -i, R n -i) = {u n ~ l ,Rn-i) is invertible. □ 



In order to apply Lemma 15.11 for the explicit monic polynomials in this paper we need to 
calculate the coefficients, which is an application of Theorem 14.51 

Lemma 5.2. Let {-R n } n >o be the monic polynomials with respect to Z on [0, 1]. Then 

™. _ y Jn F _V-^ F | y^ n(2£-j) 

n n-i 2^ 4(n + j) ^d-i 2^ 2 j ' j + 2^ A ( 2 £-j+n) j ' J+1 

j=0 v JJ j=0 i=0 v J'J 

nn _ V^ n(n-l)j(J-l) _ <A n(n - l)j 

K ^-2^32(n + j)(n + j-lf j > j - 2 ^ 8(n + j) iJ_1 

A n(n - l)(3j 2 - 6£j - 2n 2 + n - An£) 
+ ^ 16(n + j)(i - 2i - n) j ' j 



3=0 



^ n(n -I)(2Z- j) E 

2-^1 R(9P _L-n _ n\ * 



^p n{n - 1){2£ - j){2£ - j - 1) 

+ ^ 32(2£ -j+n-l){2£ + n- j) j ' J+2 

Proof. We can calculate R™-i by considering the coefficients of u n ~ 1 using the expression in 
Theorem 14.51 This gives 



( ft~\ -| \ I V L Ut i UL i Ut I (, \ / J •" J -L <Jt i Ut i <Jt lb \ / All " / 7 



using the recursive definition ( 14. 7p of [C*, [/*, V* + A°(i)] n . Note that M"_ 1 (A°(i)) is indeed 
invertible by Lemma 14.41 for irrational a. The explicit expression of the right hand side gives 
the result after a straightforward computation, since the resulting matrix is tridiagonal. 
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We can calculate R^_ 2 analogously, (-R"_ 2 )ij = 

n(n - 1) (M^Kii))- 1 ^ + n- 2)M n Q _ 1 (A«(z))" 1 {C^ + n- 1) 
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and a straightforward but tedious calculation gives the result. Note that -R™_ 2 is a five-diagonal 
matrix, since it is the product of two tridiagonal matrices. □ 

Note that even though we have used the additional degree of freedom a in the proof of 
Lemma I5.2[ the resulting expressions are indeed independent of a. 

Now we are ready to obtain the coefficients in the recurrence relation satisfied by the 
polynomials R n . 

Theorem 5.3. For any £ e |N the monic orthogonal polynomials R n satisfy the three-term 
recurrence relation 

uR n (u) = Rn+i(u) + X n R n (u) + Y n Rn-.i(u), 

where the matrices X n , Y n are given by 



2(: 



Xr, 



Y n 



-E 



t=0 



i Eij-i Ei^ f (2£ — i) E i>i+ i 



JL _i_ 



21 

E 

*=o 



4(n + z)(n + z + l) 2 A(2£ + n - i)(2£ + n - i + 1) 
n 2 (2£ + n + l) 2 



16(n + i)(n + i + l)(2£ + n - i)(2£ + n - i + I) 



Ea a. 



Proof. This is a straightforward computation using Lemma 15.11 and Lemma 15.21 The calcula- 
tion of X n is straightforward from Lemma 15.11 and Lemma 15.21 In order to calculate Y n we 
need X n R n _ 1 . A calculation shows 



X n R^_ 1 



21 

E 

j=0 



nj 2 (j - 1) E 

16(n + j-l)(n + j)(n + j + l) %3 ~ 2 



21 

E 

j=0 



nj(n + 2j + I) 



2£ 

E 

3=0 



n 
4 



Kn + j)(n + j + l) Ej ' j 1+ 

-i 2 n{2£ -i + l) 
16(n + i){n + i + 1){2£ - i + 1 + n) 

{2£-i) 2 {i + l)n 
16(21 - i + 1 + n)(2£ - % + n)(n + i + 1) 

A n(2£-z)(4£-2j+n + l) _ 
~ Z^ q^ i ~ a\(oo i ~ TTT^i.i+1 



£ 



././ 



i=o 



J(2£ + n- j)(2^ + n- j + 1)' 



2( 



+E 

i=o 



n(2£-j) 2 (2£-J + l) 



16(2£ -j + n-l)(2£-j + n)(2£ - j + n + 1) 



^iJ+2 



Now Lemma 15.11 and a computation show that Y n reduces to a tridiagonal matrix. 



□ 
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Now (14.11) and Theorem 15.31 give the three-term recurrence 

xP n {x) = P n+l {x) + (1 - 2X n )P n {x) + 4Y n P n _ 1 (x) (5.1) 

for the monic orthogonal polynomials with respect to the matrix-valued weight W on [— 1, 1]. 
The case £ = corresponds to the three-term recurrence for the monic Chebyshev polynomials 
U n . Note moreover, that hin^^oo X n = | and linij^oo Y n = i, so that the monic matrix- valued 
orthogonal polynomials fit in the Nevai class, see [TJ. Note the matrix- valued orthogonal 
polynomials P n in this paper are considered as matrix-valued analogues of the Chebyshev 
polynomials of the second kind, because of the group theoretic interpretation [T5] and Section 
but that these polynomials are not matrix-valued Chebyshev polynomials in the sense of 
13 §3]- 
Using the three-term recurrence relation (15.11) and (11. 2p we get 

4Y n # n _ 1 = f xP n (x)W(x)P n . l (xTdx 

.i (5-2) 

P n (x)W(x)(xP n ^(x)Ydx = H n 
i 

analogous to the scalar-valued case. Since Hq is determined in (II. 3p we obtain H n . 
Corollary 5.4. The squared norm matrix H n is 

(R) _ vr (n\f (21 + l)l +1 2- 2 » 

1 nh3 " l3 2(i + l)l(2£-i + l)l(n + i + l)(2l-i + n + l) 

and JH n J = H n . 

In [T51 Thm. 4.8] we have stated the three-term recurrence relation for the polynomials 
Q e n (a), a G A*, see also Section [7J of this paper. Apart from a relabeling of the orthonormal 
basis the monic polynomials corresponding to Q n are precisely the polynomials P n , see [T5J 
§6.2, (6.4)] for the precise identification 

Pd( x )n,m — ^ d Q ^d( a arccos x) -£+n,-£+m, 71, 771 G {0, 1, ... , 2£}, (5.3) 

see also Section [JJ, where T^ in (15.31) is the leading coefficient of Q d . 
Corollary 5.5. The polynomials Q e n as in [151 §4] satisfy the recurrence 

<p{a)Q E n {a) = A n Q t n+1 {a) + B n Q e n {a) + C n Q l n _ x [a) 



where 
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' (n + l)(2£ + n + 2) 

" q t^ e 2(£-q + n + l)(£ + q + n + l) 9 ' 9 ' 

r V" (l-g + l)(l + g) 

ti2(i-g + n + l)(£ + g + n + l) M_1 

(< + ? + l)(<-g) 



E 



2{£ - q + n + l){£ + q + n + 1) Em+1 



v n(2£ + n + l) 

n ^ 2(n + q + £ + l)(n - q + £ + 1) q ' r 

q t 

Proof. We use A n = T„T^ l5 B„ = T n (l - 2X„)T- 1 , C n = T„(4V;)T r ;\ and Theorem ESI 
to obtain the result from a straighforward computation. □ 

Recall from [T5j Thm. 4.8, (3.2)] that the matrix entries of the matrices A n , B n and C n 
are explicitly known as a square of a double sum with summand the product of four Clebsch- 
Gordan coefficients, hence Corollary 15.51 leads to an explicit expression for this square. 

6. The matrix- valued orthogonal polynomials related to Gegenbauer and 

Racah polynomials 

The LDU-decomposition of the weight W of Theorem 12.11 has the weight functions of the 
Gegenbauer polynomials in the diagonal T, so we can expect a link between the matrix- 
valued polynomials P n {x)L{x) and the Gegenbauer polynomials. We cannot do this via the 
orthogonality relations and the weight function, since the matrix L also depends on x. Instead 
we use an approach based on the differential operators D and E of Section [31 and because 
of the link to the matrix-valued hypergeometric differential operator as in Theorem 14.51 we 
switch to the matrix-valued orthogonal polynomials R n and x — 1 — 2u. It turns out that 
the matrix entries of P n (x)L(x) can be given as a product of a Racah polynomial times a 
Gegenbauer polynomial, see Theorem 16. 2[ 

We use the differential operators D and E of Corollary 14. 1[ and as in Section H] it is handier 
to work with the second-order differential operator D_2t — D — 2£E. By Theorem 12 .11 we have 
W(x) = L(x)T(x)L(x)\ hence Z{u) = L(l - 2u)T(l - 2u)L(l - 2m) 4 . For this reason we look 
at the differential operator conjugated by M{u) = L{\ — 2u). 

In general, for D = ^^(w) + 4-F\{u) + F (u) a second order matrix- valued differential 
operator, conjugation with the matrix-valued function M , which we assume invertible for all 
u, gives 



-, d 2 , d ( , dM- 1 \ 

M~ l DM = —M- 1 F 2 M + — M _1 FiM + 2— — F M 
du 2 du \ du J 



[M- x FoM+^—F x M+^L-FqM 

(Jj Hi Lb Lb 
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Note that differentiating M~ l M = I gives ^j^- = -M^^-M' 1 , and similarly we find 
^Mzi = -M^^M- 1 + 2M" 1 ffM- 1 ^M- 1 . We are investigating the possibility of 
M~ 1 DM being a diagonal matrix- valued differential operator. We now assume that F 2 {u) = 
w(l — u), so that M~ 1 F 2 M = u(l — u). A straightforward calculation using this assumption 
and the calculation of the derivatives of M~ l shows that M~ 1 DM = u(l — u)4-2 + J^T + T> 
with T and T\ matrix- valued functions if and only if the following equations (16. ip . (16. 2 p hold: 

FoM - d -^n-u ( l-u } d ^ = MT (6.1) 

F X M - 2u(l-u)-r- = MT t . (6.2) 

du 

Of course, T and T need not be diagonal in general, but this is the case of interest. 

Proposition 6.1. The differential operator T> = M~ l D_ 2 ^M is the diagonal differential op- 
erator 

,d 2 ( d , 
V = u 1 - u — - + — )T 1 (u + T , 



du 2 \ du 
where 

21 2£ 



T 1 {u)= l -Tl-uTl Tl = ^(2* + 3)T v , T = ^(2£ - i){2£ + i + 2)E t 



i=0 



2 

i=U 

Moreover, lZ n (u) = R n {u)M{u) satisfies 

K n V = A n {V)K n , A n {V) = A n {D)-2£A n {E). 

The proof shows that M~ 1 D a M can only be a diagonal differential operator for a = —2£. 
Note that T> is a matrix- valued differential operator as considered by Tirao, see Remark 14.21 
and [2T], and diagonality of T> implies that the matrix- valued hypergeometric 2-^i-series can 
be given explicitly in terms of (usual) hypergeometric series. In particular, we find as in the 
proof of Theorem 14.51 that 

t 



{K n {u)) k .= I 2 H t I ;uj vj 



(6.3) 



v k = (n n (o)) k ., \ n (k) = A n (v) 



kk- 



since the condition a{\Tl) <f_ — N is satisfied. 

Proof. Consider D a = D + aE, so that F 2 (u) = u(l — u) and the above considerations apply 
and F\(u) = C a — uU a , and F Q = —V a . We want to find out if we can obtain matrix- valued 
functions T\ and T satisfying (16. ip . (16. 2 p for this particular F%, F 2 and M{u) — L(l — 2u). 
Since T is diagonal, and assuming that T , T\ can be taken diagonal it is clear that taking 
the (k, l)-th entry of (16 .ip leads to 

(F ) kk M kl - ^{Tx)u ~ u{\ - u)'^- = M kl (T ) u . (6.4) 
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By Theorem 12.11 we have M k i = for I > k and for I < k 

.t / x fk\ „ fl-k,k + l + 2 

so that (16.41) has to correspond to the second order differential operator 

u{l-u)f"{u) + (c-(a + b+l)u)f'(u)-abf(u)=0, f(u) = 2 F X P^uJ (6.5) 

for the hypergeometric function. This immediately gives 

(Ti)„ = I + ^ - {21 + 3)u, (T )h - (F ) fcfe = fc 2 + 2A: - (Z 2 + 2/). 

Since (F ) fcfe = (-V a ) kk = -k 2 + (2£ + a^ff±)k - a(2£ + 2), this is only possible for a = -2£, 
and in that case 

(Ti)« = / + | - (2/ + 3)w, (T )„ = -/ 2 -21 + 2£(2£ + 2). (6.6) 

It remains to check that for a = —21 the condition (16. 2p is valid with the explicit values 
(16. 6p . For a = —2£ the matrix- valued function F± is lower triangular instead of tridiagonal, 
so that (16.21) is an identity in the subalgebra of lower triangular matrices. With the explicit 
expression for M we have to check that 

({- + k)-u{3 + 2k))M k i-kM k - lf i - 2u{l-u)-r^ = M kl (- + l-u{3 + 2l)) 
which can be identified with the identity 

dC (l+l) 
(1 - x 2 ) -^-{x) = (k + l + l)Cti\(x) - x(k - l)ctl\x). 

In turn, this identity can be easily obtained from [Tj (22.7.21)] or from [TTj (4.5.3), (4.5.7)]. □ 

Since R n and M are polynomial, Proposition 16.11 and the explicit expression for the eigen- 
value matrix in Corollary 14.11 imply that (TZ n ) k j is a polynomial solution to 

u(l - u)f"(u) + ((j + |) - u(2j + 3)) f(u) + {21 - j)(2£ + j + 2)/(i») 



-n(n - 1) - n(2£ + 3) + A;(2£ - jfc) + 2n(£ - k) - (2£ + 2)(k- 2£)\ f(u) 
which can be rewritten as 

«(1 - «)/"(«) + ((j + ^) - «(2j + 3)) /'(«) - (j - k - n)(n + k + j + 2)f(u) = 

which is the hypergeometric differential operator for which the polynomial solutions are 
uniquely determined up to a constant. This immediately gives 

t? i \ r ^ r? fJ-k-n,n + k + j + 2\ 

T<-n{u)kj = c kj {n) 2-Fi I 3 ;u\ . (6.7) 

for j — k — n < and TZ n (u) k j = otherwise. The case n = corresponds to Theorem 12.11 
and we obtain c k j(0) = ( .). It remains to determine the constants c k j(n) in (16.71) . 
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First, switching to the variable x, we find 

(V n {x)) k . = (P n {x)L{x)) k . = (-2)»qy(n) Qj + 2ul\ C ^-i {x) (6 ' 8) 

so that by ( 16 .8p the orthogonality relations (jl.2p and (I2.2p give 

2M(n+fc) 

5 nm (H n ) kl = (-2) n+m J2 %-( fl )%H c i( l ) 

j=0 

(n + k-j)\ v^r(j + |) 

X (2j + 2) n+k _ J dn+k > m+l (n + k + l)j\- 

Using the explicit value for Cj (£) as in Theorem 12.11 and Corollary 15.41 we find orthogonality 
relations for the coefficients c/~j(n): 

2M(n+fe) 
-In - *~~^ 



x H n )kk^ " n 5 nm — 2_j c kj{n)Ck+ 



m—n,j 



m 



j=° (6.9) 

(j!) 2 (2j + l)(2l + j + l)!(2l-j)!(n + fc-j)! 
(n + k + j + 1)! (n + k + 1) (2£)\ 2 

Note that we can also obtain recurrence relations for the coefficients Ckj(n) using the three- 
term recurrence relation of Theorem 15.31 

Theorem 6.2. The polynomials lZ n (u) = R n (u)M(u) satisfy 
v ( s ( v ^ (-2i)j(-k-n)j f-j,j + l,-k,-2£-n-l 

^ fj-k-n, n + k + j + 2 

X2 M ;+i 



n j— n 



n\ {21 + 2) n 



mt/i H n (u)kj = /or j — k — n > and 

<*.<») = <-D"4 (fc + !).(«_* + !). 

We view Theorem 16.21 is an extension of Theorem 12.11 but Theorem 12.11 is instrumental in 
the proof of Theorem 16.21 Since the inverse of M(u), or of L(x), does not seem to have a 
nice explicit expression we do not obtain an interesting expression for the matrix elements 
of the matrix-valued monic orthogonal polynomials R n (u) or of P n (x). Note also that the 
case £ = gives back the hypergeometric representation of the Chebyshev polynomials of the 
second kind U n , see (14. 9p . 

Comparing with (12. 5 p we see that the 4.F3-series in Theorem 16.21 can be viewed as a Racah 
polynomial Rk(X(j); —2£ — 1, —k — n — 1, 0, 0), respectively _Rfc+n-m(A(j); — 2£ — 1, —k — n — 
1,0,0), see (12. 5p . where the N of the Racah polynomials equals 2£ in case 2£ < k + n and 
N equals k + n in case 2£ > k + n. Using the first part of Theorem 16.21 we see that the 



MATRIX- VALUED ORTHOGONAL POLYNOMIALS II 21 

orthogonality relations ( 16. 9 p lead to 

(H) 2-^S -1 2i+1 \ r M |2^V fc) (2j + l)(-2^(-n-fc) 3 - 

1 n)kk nm ~ 2(n + k + iy l k0[ )l f^ Q (2£ + 2) j (n + k + 2) j (6.10) 

x R k (X(j); -2£ - 1, -k - n - 1, 0, 0) it! fe+n _ m (A(j); -2£ - 1, -fc - n - 1, 0, 0), 

which corresponds to the orthogonality relations for the corresponding Racah polynomials, 
see [21 p. 344], [HJ §1.2]. From this we find that the sum in (I6.10p equals 

(2£ + l)(n + k + 1) 



Hence, 



Snm (2£+\+n-k) 



, ( s { 2 (H) 2n 2(n + k + l)(2£+l + n-k) 
\Cko{n)\ — {H n ) kk 2 



7i (2^ + 1; 



2 



A _ 2n (n\Y(2£ + 2)l 



(6.11) 



(k + l)l(2£-k + l)l 

using Corollary 15.41 

We end this section with the proof of Theorem 16. 2[ The idea of the proof is to obtain a 
three-term recurrence for the coefficients Ckj{n) with explicit initial conditions, and to compare 
the resulting three-term recurrence with well-known recurrences for Racah polynomials, see 
[2], [H], [E]. The three-term recurrence relation is obtained using the first-order differential 
operator E and the fact that the lZ n , being analytic eigenfunctions to T>, are completely 
determined by the value at 0, see Remark 14.21 

Proof of Theorem \6.B. Since the matrix- valued differential operators D and E commute and 
have the matrix- valued orthogonal polynomials R n as eigenfunctions by Corollary 14. 1\ we see 

that £ = M~ l EM satisfies 

£K n = A n (£)n n , A n {£) = A n (E), £V = V£ (6.12) 



Moreover, in the same spirit as the proof of Proposition 16.11 we obtain 

S x {u) = n{l-u)^ w - w - TY) E^ 1 - y^—g-E^ (6 . 13) 



i=0 v /\ / j =0 

21 .o /„/,... -s 2£ 



S (u) = (1 - 2u) ^ _ E^ + ^ 

i=0 v ' i=0 

by a straightforward calculation. 

Define the vector space of (row-)vector valued functions 

V(\) = {E analytic at u = | FV = XE}, 



E; 
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and v. V(X) — > C 2e+l , F i— y F(0), is an isomorphism, see Remark 14.21 and [21J. Because of 
(I6.12p we have the following commutative diagram 

V{\) — *—> V{\) 



£21+1 N ( X ) ) (Q2£+l 

with N(X) a linear map. In order to determine N(X) we note that F e V(A) can be written 
as, cf (16. 3p . 






F,( M )= 2 ffJ *'" u ;u F(0) 



2-I ' /j 



dFj (c\\ i?ic\\t\ ^ ^^('l ^ ^l 



N{X) = (X-T )(^Tl)- 1 S 1 (0) + S (0) 



, 2( 



so that -^f-(O) = F(0)(A — T )(|T 1 1 ) x by construction of the 2 -Hi-series, see Remark |4"72"1 Now 
(IBTT21) gives 

acting from the right on row- vectors from C +1 . 

By Proposition 16. II we have that the /c-th row ((7£ n ) &;(•)) _ of R n is contained in V(X n (k)), 
see ( 16. 3p . On the other hand, the k-th row of lZ n is an eigenfunction of £ for the eigenvalue 
jUn(fc) = A n (£) kk . Since ^f (("^n)fej) i=0 J = (c fei (n))^£ we see that the row-vector c fc = 
(cfcj(^))jlo satisfies CkN(X n (k)) = fi n (k) Cj., which gives the recurrence relation 

(i + k + n+l)(i-k-n-l)(2£-i + l) 



(2i + 1) 



-c fc>i _i(n) + (z(i + 1) - 4£(£ + l))c fcii (n) 



+ ( ' + 1) ( ^y + 2) c M+1 (n) = (-2n(l - fc) + (2£ + 2)(fc - 2*)) Cfcji (n), (6.14) 

with the convention c&_i(n) = 0. Note that Cjfe(O) = ( ■) indeed satisfies ( I6.14p . Comparing 
( I6.14p with the three-term recurrence relation for the Racah polynomials or the corresponding 
contiguous relation for balanced 4-F 3 -series, see e.g. [2, p. 344], [HI §1.2], gives 

, n / v i^ (-^)j(- fe - w )j F /-J,J + l,-fc,-2£-n-l 

and Ckjiji) = for j > k + n. 

It remains to determine the constants c/~o(n), and we have already determined their absolute 
values in (16. lip by matching it to the orthogonality relations for Racah polynomials. From 
the three-term recurrence relation Theorem 15.31 we see that the constants Ckj{n) are all real, 
so it remains to determine the sign of Cko (n) . Theorem 15.31 gives a three-term recurrence for 
lZ n (u), and taking the (k, 0)-th matrix entry gives a polynomial identity in u using (16.71) . Next 
taking the leading coefficient gives the recursion 

(n + k + 2) {21 -kf 

Cfco(n + 1} = - 4(n + fc + l) Ck ° {n) + A(2i + n-k)(2i + n-k + l) ° k+1 ' M 
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and plugging in Ckoiji) = sgn.(cko(n))\cko(n)\ and using the explicit value for |cfco(n)| gives 

sgn(c fc0 (n + 1)) (n + 1){2£ + n + 2) = 

- sga(c k0 {n)){n + k + 2) (2/ - k + n + 1) + sgn(c fc+ i, (n))(2£ -k){k + l). 

This gives sgn(cfc (^)) = sgn(cfc + i )0 (n)) for the right hand side to factorise as in the left hand 
side, and then sga(cko(n + 1)) = — sgn(cfc (n)). Since Cfc (0) = 1, we find sgn(cfc (n)) = 
(-1)". □ 



Remark 6.3. Theorem 16.21 can now be plugged into the three-term recurrence relation for 1Z n 
of Theorem I5.3[ and this then gives a intricate three-term recurrence relation for Gegenbauer 
polynomials involving coefficients which consist of sums of two Racah polynomials. We leave 
this to the interested reader. 

Remark 6.4. We sketch another approach to the proof of the value of Cko{n) by calculating 
the value Ck t 2e{ n ) i n case k + n > 2£ or Ck t n+k{ n ) i n case k + n < 2£. For instance, in case 
k + n > 2£ we have 

(H n (u)) kj2 e = {R n {u)M(u)) k ,2i = (R n (u))k,2e = (Rn(u)) 2 e-k,o 

using that M is a unipotent lower-triangular matrix-valued polynomial and the symmetry 
JR n {u)J = R n (u), see [H2 §5]. Now the leading coefficient of the right hand side can be 
calculated using Theorem 14.5} and combining with (16. 7p . the value Ck,2e{ n ) follows. Then the 
recurrence (I6.14p can be used to find Cko{n). 

7. Group theoretic interpretation 

The purpose of this section is to give a group theoretic derivation of Theorem 13.11 comple- 
menting the analytic derivation of [H3, §7]. For this we need to recall some of the results of 



7.1. Group theoretic setting of the matrix- valued orthogonal polynomials. In this 
subsection we recall the construction of the matrix- valued orthogonal polynomials and the cor- 
responding weight starting from the pair (SU(2) x SU(2), SU(2)) and an SU(2)-representation 
T . Then we discuss how the differential operators come into play and what their relation is 
with the matrix-valued orthogonal polynomials. The goal of this section is to provide a map 
of the relevant differential operators in the group setting to the relevant differential operators 
for the matrix- valued orthogonal polynomials in Theorem 17.81 

Let U = SU(2) x SU(2) and K = SU(2) diagonally embedded in U. Note that K is the set 
of fixed points of the involution 9 : U —$■ U : (x,y) H- (y,x). The irreducible representations 
of U and K are denoted by T f ~ 1,£2 and T e as is explained in [TSJ §2]. The representation space 
of T e is denoted by H which is a 2£+ 1-dimensional vector space. If T e occurs in T £l,£2 upon 
restriction to K we defined the spherical function $^ l2 in [151 Def. 2.2] as the T £ -isotypical 
part of the matrix T £l ' £2 . Let A C U be the subgroup 

A ={ at =(( e T e -/0'( e T e ^))'°^ t<47r 
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and let M = Zk(A). Recall the decomposition U = KAK, [12, Thm. 7.38]. The restricted 
spherical functions $1 ^\a take values in EiadM(H e ), see [15, Prop. 2.3]. Since EndM(-f^ £ ) — 
C 2e+l this allows allows us to view the restricted spherical functions as being C 2£+1 -valued. 
The parametrization of the [/-representations that contain T e indicates how to gather the 
restricted spherical functions. Following [151 Fig- 3] we write (^1,^2) — C(d,h) with ((d,h) = 
(±(d + £ + h),±(d + £-h)). Here d e N and h e {-£, -£ + l,...,£}. We recall the definition 
of the full spherical functions of type £, [T51 Def. 4.2]. 

Definition 7.1. The full spherical function of type £ and degree d is the matrix-valued function 
§ d : A — y End(C 2 ^ +1 ) whose j-th row is the restricted spherical function $1 £a with [£^£2) = 
C(d,j). 

The full spherical function of degree zero has the remarkable property of being invertible 
on the subset A reg := {a t \t G (0, n) U (ir, 2n) U (2ir, 3n) U (3tt, 47r)}, which was first proved by 
Koornwinder [TOJ Prop. 3.2]. The invertibility follows also from Corollary I2.3I Let = $°/ 2 l , 2 
be the minimal non-trivial zonal spherical function [T5"| §3]. Together with the recurrence 
relations for the full spherical functions with <fi [T5l Prop. 3.1] this gives rise to the full spherical 
polynomials [151 Def. 4.4]. 

Definition 7.2. The full spherical polynomial Q e d : A —y End(C 2£+1 ) is defined by Q e d {a) = 
(H(a))-^ d (a)- 

The name full spherical polynomial comes from the fact that the Q d are polynomials in 0. 
The full spherical polynomials Q d are orthogonal with respect to 

(Q, P) v t = f Q(a)V e (a)P(a)da, V e (a t ) = ($^(a 4 ))*$^(a 4 ) sin 2 1, 

J A 

see [13 Cor. 5.7]. 

In [13 §5] we studied the weight functions V extensively. It turns out that the matrix 
entries are polynomials in the function 0, apart from the common factor sint. Upon changing 
the variable x = 0(a) we obtain the following system of matrix- valued orthogonal polynomials. 

Definition 7.3. Let R d : [0, 1] -)> End(C 2£+1 ) be the polynomial defined by R d ((f>(a)) = Q e d (a). 
The degree of R d is d. The polynomials are orthogonal with respect to 



(R,P) w t= J R(x)W e (x)P(x)dx, 



where W e is defined by W e (<ft(a))d<f) = V e (a)da. 

The weight W e (x) from Definition 17.31 is the same as the weight defined in (11.11) where we 
have to bear in mind that the basis is parametrised differently. The matrix- valued polynomials 
R d correspond to the family {Pd}d>o from Theorem 13.11 by means of making the R d monic. 
Given a system of matrix-valued orthogonal polynomials as in Definition 17.31 it is of great 
interest to see whether there are interesting differential operators. More precisely we define 
the algebra V(W e ) as the algebra of differential operators that are self-adjoint with respect 
to the weight W l and that have the R d as eigenfunctions. We define a map that associates to 
a certain left invariant differential operator on the group U an element in T>(W). 
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Before we go into the construction we observe that the spherical functions may also be 
defined on the complexification A c , using Weyl's unitary trick. Indeed, all the representations 
that we consider are finite dimensional and unitary, so they give holomorphic representations 
of the complexifications U c and K . 

A great part of the constructions that we are about to consider follows Casselman and 
Milicic [5], where the differential operators act from the left. In this section we follow this 
convention, except that we transpose the results at the end in order to obtain the proof of 
Theorem 13.11 where the differential operators act from the right. 

Let U(u c ) be the universal enveloping algebra for the complexification u of the Lie algebra 
u of the group U = SU(2) x SU(2). Let 9: U(u c ) ->■ U(u c ) be the flip on simple tensors 
extending the Cartan involution 9: su(2) x su(2) -» su(2) x su(2), (X,Y) (->■ (Y,X). Recall 
t = su(2) is the fixed-point set of 9. Let £/(u c )* denote the subalgebra of elements that 
commute with t c . Let 3 denote the center of U(i c ). 

Lemma 7.4. U(u c f C = 3 ® 3 ® 3- 

Proof. From [13, Satz 2.1 and Satz 2.3] it follows that f/(u c ) {C = 3 <8>z(j) (3 ® 3) where j C u c 
is the largest ideal of u c contained in £ c . Since j = the result follows. □ 

Proposition 7.5. The elements of the algebra U(u c Y have the spherical functions $| i2 as 
eigenf unctions. This remains true when we extend $| £a to U c . 

Proof. See [221 Thm. 6.1.2.3]. The second statement follows from Weyl's unitary trick. □ 

The spherical functions $| £a have T^-transformation behaviour: 

¥ ll/2 {k lU k 2 ) = T^^lfeHT^fe) (7.1) 

for all ki,k 2 G K and u G U, see [T5J Def. 2.2]. Let C(A) denote the set of continuous 
(C-valued) functions on A. Casselman and Milicic [Sj define the map 

IT*: U{u c f ->• C{A) ® U{a c ) ® End(End M (i/ £ )) 
and prove the following properties [5] Thm. 3.1, Thm. 3.3]. 

Theorem 7.6. Let F : U — > End(/7 £ ) be a smooth function that satisfies (j7.ip . T/ien (_D.F)|yi = 
II^(£))(i r '|A) /o?" ffl^ D G f/(u c ) t . Moreover, Hi is an algebra homomorphism. 

We call ILj(£)) the T £ -radial part of D G t/(u c )* . In particular we have 

n e (D)($l iji2 \ A ) = \ e D/l/2 $i >e2 \A: Ai, AA G c. 

Upon identifying End a/ (H e ) = C 2e+l we observe that we may view ILj(£)) as differential 
operator of the End(C 2£+1 )-valued functions that act on from the left. In particular, let 
C(A, End(C 2 ^ +1 ),T £ ) denote the vector space generated by the $^, d > 0. The following 
lemma follows immediately from the construction. 

Lemma 7.7. Let D G £/(u c )* be self-adjoint and consider U-t(D) as a differential operator 
acting on C(A, End(C 2£+1 ),T^) from the left. Then He(D) is self- adjoint for (•, -) v e. 

Definition 7.8. Let f : U(u ) — > V(W e ) be defined by sending D to the conjugation of the 
differential operator Tl^D) by $q followed by changing the variable x = <p(a). 
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The map /: £/(u c ) { — > T>(W ) is an algebra homomorphism. It gives an abstract descrip- 
tion of a part of T>(W ). Note that / is not surjective because in [T5], Prop. 8.1] we have found 
a differential operator that does not commute with some of the other. However, the algebra 
U(u c y is commutative by Lemma [7.41 

By means of Lemma 17.41 we know that f/(u c ) t has Qi = ty <8> 1 and Q2 = 1 <8> ty among 
its generators, where f2f e 3 is the Casimir operator. In the following subsection we calculate 
/(Qi+f^) and f(Qi— f^) explicitly. Upon transposing and taking suitable linear combinations 
we find the differential operators D and E from Theorem 13.11 



7.2. Calculation of the Casimir operators. The goal of this subsection is to calculate 
/(fi) and f(fi?) where / is the map described in Definition 17.81 and where Q = Qi + f2 2 and 
Q' = Qi — Q 2 - We proceed in a series of six steps. (1) First we provide expressions for 
the Casimir operators Q and Q' which (2) we rewrite according to the infinitesimal Cartan 
decomposition defined by Casselman and Milicic [5, §2]. These calculations are similar to those 
in [221 Prop. 9.1.2.11]. (3) From this expression we can easily calculate the T £ -radial parts, see 
Theorem 17.61 The radial parts are differential operators for Endjvf(# )-valued functions on 
A. At this point we see that we can extend matters to the complexification A c of A as in (5J 
Ex. 3.7]. (4) We identify EndM(H e ) — C 2£+1 and rewrite the radial parts of step 3 accordingly. 
(5) We conjugate these differential operators with $0 and (6) we make a change of variables to 
obtain two matrix- valued differential operators f(£i) and /(O'). Along the way we keep track 
of the differential equations for the spherical functions. Finally we give expressions for the 
eigenvalues A d and T d of f(fl) and f(fi') such that the full spherical polynomials Q d are the 
corresponding eigenfunctions. Following Casselman and Milicic [5-, §2] the roots are considered 
as characters, hence written multiplicatively. 

(1). First we concentrate on one factor K = SU(2), with Lie algebra t and standard Cartan 
subalgebra t. The complexifications are denoted by £ c , t c and we use the standard basis 

H - ( l M F -f° A F if 00 

~ [0 -1) ' a 2 V° °y ' 2 V 10 

for ^ c . The Casimir of K is given by fl? = \H 2 + 4 {E a E a ~i + E a ~iE a }. It is well-known that 
the matrix-elements of the irreducible unitary representation T l of SU(2) are eigenfunctions 
of the Casimir operator fi t for the eigenvalue |(£ 2 + £), see e.g. [121 Thm. 5.28]. The roots of 
the pair (t c ®t € , t c ©t c ) are given by R = {(a, 1), (a -1 , 1), (l,a), (l,a -1 )}. The positive roots 
are choosen as R + = {(a, 1), (1, a -1 )}, so that the two positive roots restrict to the same root 
R(u c , a € ) which we declare positive. The corresponding root vectors are Ei a ^ = (E a , 0), etc. 
Define 

E = (E a ,0)(E a -i,0) + (E a -i,0)(E a ,0), 

Then we have Qi = \(H, 0) 2 +4E and f2 2 = 0(Qi). In particular, the spherical function $1 i2 
is an eigenf unction of Qi for the eigenvalue |(£? + IA for i = 1, 2. 
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We have (H, 0) = §((#, -H) + (if, /J)) and (0, H) = \{{H,H)- (H, -H)) and from this 
we find in U(u c ) 

n = n, + n 2 = 1 -(h,h) 2 + -(h,-h) 2 + a(e + 6(e)), 

1 4 (7-2) 

Q' = Qx-n 2 = ~(H,-H)(H,H)+4(E-6(E)). 

(2). Following Casselman and Milicic [5J §2] we can express Q and Q' according to the 
infinitesimal Cartan decomposition of U(u c ). Let G R and denote X^ = Ep + 9{Ep) G t . 
Denote y a = Ad(a _1 )y for a G A In [SJ Lemma 2.2] it is proved that the equality 

(l-/3(a) 2 )X /3 = /3(a)(i^-/3(a)£ /3 ) 

holds for all a G A reg . This is the key identity in a straightforward but tedious calculation to 
prove the following proposition which we leave to the reader. 

Proposition 7.9. Let a G A reg and (3 G R + . Then 

Q = re «* - H ^ 2 + ^ ^ 2 ) - m^mr {xm - 1+ 



X a ^X% + XpXp-i + Xp-iX p - (/3(a) + (5{a)- l ){X*X^ + X^X fi )} 

l/3(a) + /3(a 
4/3(a)-/3(a 






and 



fi - ~(H, H)(H, -H) + j— _ /3(a) _ 1 (if, if) 

+ fta) - ft(a)-i {X e- lXfi ~ X ^" l} ' (7 ' 4) 

The calculation of (17. 3p is completely analogous to that of [22| Prop. 9.1.2.11] and it is clear 
that (17.31) is invariant for interchanging (3 and /3 _1 . The expression in (17.41) is also invariant for 
interchanging j3 and /3 _1 albeit that it is less clear in this case. In either case the expressions 
(17. 3p and (I7.4p do not depend on the choice of j5 G f? + . 

(3). Following Casselman and Milicic [3 §3] we calculate the T £ -radial parts of Q and 
Q'. This is a matter of applying the map LT^ from Theorem 17.61 to the expressions (17. 3p and 
(I7.4p . At the same time we note that the coefficients in (17. 3p and (17.41) are analytic functions 
on A reg . They extend to meromorphic functions on the complexification A c of A which we 
identify with C x using the map 

w \ / w^ 1 



a : C x — > A : w i— > a(w) — i i ,, _ i . i 

w J y u> 

Under this isomorphism the differential operator (ii, —if) translates to Wt-. To see this let 
g : A c — » C be holomorphic and consider (ff, —H)g(a(w)) which is equal to 

(H, -H)g(a(w)) = I — #(a(e*w)) \ = w—(g o a)(w). 
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Following [5], [22] we find the following expressions for the T^-radial parts of Q and Q'; 
^ /^\ 1 ( d\ lw 2 + w~ 2 d 1 m » /TT , 9 

{tHEJTHE^) . +T'(E a -,)T'(E a ) . + 



(w 2 — w 2 ) 2 



>T\E a )T%E a -,) + •T e (E a - 1 )T e (E a )} + 

{T\E a ) . T<(£ a -i) + T\E a -,) . T £ (E Q )} , (7.5) 



w 2 + w 2 
(w 2 — w~ 2 ^ 2 



and 



d , lw 2 + w~ ., 



TWO) = ~T\H)wf- + j 2 _ 2 T\H)+ 
8 dw 4w z — w 2 



w 2 — w~~ 2 



{T\E a -,).T\E a )+T\E a ).T\E a -,)} (7.6) 



where the bullet (•) indicates where to put the restricted spherical function. The matrices 
T e (E a ) and T e (H) are easily calculated in the basis of weight vectors. Note that T (E a -i) = 
JT e (E a )J. We give the entries of T e (E a ) in the proof of Lemma [7.111 

The following proposition is a direct consequence of Theorem 17.61 and Proposition 17.51 

Proposition 7.10. The restricted spherical functions are eigenf unctions of the radial parts 
of Q and fi' , 

n £ (fi)($?'V) = ^i + 4 + ^ + 4)<&?''V, 

(4). The spherical functions Q\ £z restricted to the torus A c take their values in End.M(H ) 
and this is a 11 + 1-dimensional vector space. We identify 

End M (#0 -> C 2m :Y^Y up 

to obtain functions ($1 f J^c) up . The reason for putting the diagonals up is that we want 
to write the differential operators as differential operators with coefficients in the function 
algebra on A with values in End(C 2 ^ +1 ) instead of the way ILj(f2) and ILj(f2') are defined. The 
differential operators that are conjugated to act on C 2f+1 -valued functions are also denoted 
by (") up - The differential operators ( 17. 5p and (17. 6p that are defined for EndM(-£^)-valued 
functions conjugate to differential operators Ue(Q) up and Ue(Q') np for C 2 ^ +1 -valued functions. 
All the terms except for the last ones in (17. 5p and (17. 6p transform straightforwardly. 

Lemma 7.11. The linear isomorphism End.M{H ) — ¥ C 2e+1 : D \— > D up conjugates the linear 
map End M (H e ) ->■ End M (#*) : D i-» T e {Ej)DT e (E a -x) to C 2e+1 ->■ C 2e+1 : D np h-> C e D up , 
where C G End(C 2 ^ +1 ) is the matrix given by 

c £ p , = \(i + j)(i - j + i)W, t< Pi j<e. 
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Likewise, D i— )■ T e (E a -i)DT e (E a ) transforms to D up i-> JC f JD up , where J is the anti- 
diagonal defined by Jij = 8^-j with —C. < i,j < t. 

Proof. Working with the normalized weight-basis as in [161 §1] we see that T e (E a ) is the 
matrix given by 

T*(v \ A i + l + 1 / (^-^-2)!(f + z"+2J! 
[J^ahj o M+ i 2 y (£_i_ l)!(£ + i + l)! 

and T e (E a -i) = JT e (E a )J. The lemma follows from elementary manipulations. □ 

We collect the expressions for the conjugation of the differential operators (17. 5p and (17.61) 
by the linear map V i-> V up where we have used Lemma 17.111 

[ {T i (E a )T\E a -i)+T\E a -i)T\E a )} + 



(w 2 — w 



-2\2 



W 2 + W 2 
(w 2 — W~ 2 ) 2 



{JC'J + C 1 }, (7.7) 



8 aw 4w 2 — w 2 



+ 2 2 2 {JCV-^}. (7.J 



a 



The differential operators ( 17. 7p and ( 17. 8p also act on the full spherical functions $ d ' . Col- 
lecting the eigenvalues of the columns in $ d ' in diagonal matrices we obtain the following 
differential equations: 

n £ (^i + fi 2 ) up ^ = svU (7.9) 

n,(fli - fi 2 ) up $d = $ d r d , (7.io) 

where (A d ) pi = ^(d 2 + j 2 + 2d{£ + 1) + £(£ + 2)) and (T d ) pj = \8 Ptj j{i + d+l). For further 
reference we write 

n,(fti + aO up = a 2 (w)^ + ai (w)^- + a (w), (7.11) 

n £ (fi!-fi 2 ) Up = 6 1 ( U ;)A + 6o ( u; ). (7.12) 

dw 

(5). Recall from Definition 17. 1 1 that the full spherical polynomials Q d ' are obtained from the 
full spherical functions $ d ' by the description Q d ' = ($ ' ) _1< ^d • We conjugate the differential 
operators (17. 7p and (17. 8p with $ to obtain differential operators to which the polynomials Qd 
are eigenfunctions. We need a technical lemma. 
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Lemma 7.12. Let a e : C x -)■ End(C 2£+1 ) be defined by a e (w) = £(w 2 + uj- 2 )I + S £ where S e 
p,j 



is defined by (S £ ) = -(£ - j)S p - jA - (£ + j)8j- p ,i. Then 



\w{w 2 - w- 2 )^-^\w) = ^ t (w)a\w). (7.13) 

Let v* : C x -). End(C^ +1 ) be defined by v e (w) = |^ (^^ iag + u£), where ([//J. . = 
{-21 + 2j)S i>j+1 + {2£ + 2j)S i+ld and (U^ g ) . . = -2z^-. Tften 

b 1 {w)^\a{w)) = ^(«(»))^H- (7-14) 

Proof. The matrix coefficients of $ ' (a(w)) are given by 

mm(t~p,i-j) w 4r-2e+2p+2j 

£-*< r \(£-p- r )\(£- j — r)!(p + ?' + r)!' 

r=max(0,— p— j) 

see [TB| Prop. 3.2]. The matrix- valued function bi(w) is equal to the constant matrix ^T e where 
T e (H)ij = 2Sijj. We can now express the matrix coefficients of the matrices in (17.13P and 
( I7.14p in Laurent polynomials in the variable w and comparing coefficients of these polynomials 
shows that the equalities hold. □ 

Definition 7.13. Define tt e = (^o'*)" 1 ° IL(fi) up o $J* and A e = (^q'*)" 1 ° IL(ft') up ° $o*. 
Theorem 7.14. The differential operators fli and A^ are given by 

A, = A">)^ + To (7.17) 

Proof. This is a straightforward calculation using the expressions ( 17. lip and (j7.12p . bearing 
in mind that the coefficients are matrix-valued. In both calculations the difficult parts are 
taken care of by Lemma 17.121 □ 

i i 
(6). The elementary zonal spherical function $q' 2 is denoted by and we have <p{a{w)) = 

\{w 2 + w~ 2 ). In this final step we note that the differential operators f^ and A^ are in- 
variant under the maps w i— > —w and w \-¥ w~ l . This shows that the differential opera- 
tors can be pushed forward by <fi o a to obtain differential operators on C in a coordinate 
z = <f)(a(w)). Using the identities w-j^{h o <p){w) = {w 2 — w~ 2 )h'{cf){w)), {w-£^) 2 {h o <f>)(w) = 
{w 2 - w- 2 ) 2 h"((f)(a(w))) + 2{w 2 + w W2 )h'((f)(a(w))) and {w 2 - w~ 2 ) 2 = 4{(j){a{w)) 2 - 1) we 
transform (I7.16P and (17.171) into 

A, = I (2'UL* + K,) j; + To (7.19) 
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Recall that the End(C 2 ^ +1 )-valued polynomials R d ' are defined by pushing forward the 
End(C 2 ^ +1 )-valued functions Q d ' over (ft o a, see Definition 17.31 

Theorem 7.15. The members of the family {Rj }d>o of Eiad(C 2e+1 ) -valued polynomials of 
degree d are eigenf unctions of the differential operators Q( and Ag with eigenvalues A^ and Yd 
respecively. The transposed differential operators (QeY and (A^) f satisfy 

-4(n e y + 2(£ 2 + £) = D, 
-l{A t ) t -(l + l) = E, 



where D and E are defined in Theorem \3.1[ 

Proof. The only things that need proofs are the equalities of the differential operators. These 
follow easily upon comparing coefficients where one has to bear in mind the different labeling 
of the matrices involved in the two cases. □ 

Note that the differential operators D and Q( are invariant under conjugation by the matrix 
J, where Jij = 8^-j. The differential operator A^ is anti-invariant for this conjugation. The 
differential operator E does not have this nice property. 

Appendix A. Proof of Theorem 12.11 

The purpose of this appendix is to prove the LDU-decomposition of Theorem 12. II We prove 
instead the equivalent Proposition 12.21 an d we start with proving Lemma 12.71 

Note that the integral in Lemma I27T1 is zero by (12. 2 \ in case t > m, since C m _ k (x)U n+m -2t{x) 
is a polynomial of degree n + 2m — k — 2t <n — k. 

We start by proving Lemma 12.71 in the remaining case for which we use the following well- 
known formulas for connection and linearisation formulas of Gegenbauer polynomials, see e.g. 
|2l Thm. 6.8.2], pU Thm. 9.2.1]; 



Ln/2j 



rW , v v- (7 ~ P)k(l)n-k ( {3 + n-2k \ m 
° n {X) ~ ^ k\ (/3 + l) B _ fc I ) C - 2fcl 



fe=0 



;r 



rt<*) f „\r>(«)(^ _sr^ ( n + m - 2k + a )( n + m - 2fc)!(a) fc (A.l) 

C n [x)C m [x) - ^ ________ 

(«)n-jt(a)m-fc(2a) n+m -fc (q) , . 

x ( n -k)\(m-k)\(a) n+m . k (2a) n+m . 2k "+™~ 2feW 



Proof of Lemma\2.1\ We indicate the proof of Lemma \2.1\ so that the reader can easily fill in 



the details. Calculating the product of two Gegenbauer polynomials as a sum using the lineari- 
sation formula of (1A.1|) and expanding the Chebyshev polynomial C/ n+m _2t(x) = C, 



(i) 

n+m— 2£ 



X 



in terms of Chebyshev polynomials with parameter k + 1 using the linearisation formula of 
(lA.ip , we can rewrite the integral as a double sum with an integral of Chebyshev polynomials 
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that can be evaluated using the orthogonality relations (12.21) reducing the integral of Lemma 
12.71 to the single sum 

min(t,m— k) . . ., . ,, . , . 

^ (m + n - k + J - 2r) (k + l) r (k + l) n _ k _ r {k + l) m -k-r 



(m + n — k + 1 — r) r\ (m — k — r)\ (n— k — r 
(2/c + 2) m+n _ 2fc -r(-^)fc+r-t(^ + m-t-k-r)\ y/7tT(k + 2 

X 



3> 



(k + l) m+n _ 2fe _ r (A; - t + r)\ (k + 2)„ +m _ t „ fc _ r (k + 1) T(k + 1) ' 

min(i 
r=0 



Assuming for the moment that k > t, so the sum is X^o m • Then this sum can be written 



as a very-well-poised 7-F 6 -series 

^T(k + |) (k + l) m _ fc (k + l) n _ fc (2fc + 2) m+n _ 2fc (-fc) fc _ f (n + m - t - fc)! 



(A; + 1) T(k + 1) (m-k)\ (n - k)\ (k + l) m+n „ 2k (k - t)\ (k + 2) n+m ^ k 
\{k — m — n+ l),k + l,k — m, k — n, k — m — n — 1, —t,t — m — n — 1 



|(/c — m — n — 1), — m, — n, —m — n — l,k — t + 1, —n — m + k + t 

Using Whipple's transformation j2j Thm. 3.4.4], [U §4.3] of a very- well-poised 7_F 6 -series to a 
balanced 4F 3 -series, we find that the 7.F 6 -series can be written as 

(k — m — n)t {—t)t j-, (— k, k+ 1, — t, t — m — n — 1 

4-F3 



(A; — £ + l)t (— m — n — 1)* \ — ra, — m, 1 

Simplifying the shifted factorials and recalling the definition of the Racah polynomials (12. 5p 
in terms of a balanced 4.F3-series gives the result in case k > t. 

In case k < t we have to relabel the sum, which turns out again to be a very-well-poised 
7_F6-series which can be transformed to a balanced 4_F3-series. The resulting balanced 4F3- 
series is not a Racah polynomial as in the statement of Lemma 12. 7\ but it can be transformed 
to a Racah polynomial using Whipple's transformation for balanced 4_F 3 -series |2] Thm. 3.3.3]. 
Keeping track of the constants proves Lemma [2 .71 in this case. □ 

As remarked in Section |2j Theorem 12.11 follows from Proposition 12.21 In order to prove 
Proposition 12. 21 we assume act(m, n) to be known ( II. ip and to find (3k{m, n). Given the explicit 
a t (m,n), and multiplying by \/l — x 2 C/ n+m _ 2 t(x) and integrating we find from Lemma [2T7I 

m 

a t (m,n)- = )^ k (m, n)C k (m, n)R k (X(t); 0, 0, -m - 1,-n- 1) (A. 2) 

k=0 



where 



C k (m,n) 



^T(k +§)(* + l) m _ fc (k + l) n _ k {-If (2k + 2) m+n _ 2k (k + 1)! 



(A; + l) (m-k)\ (n-k)\ (n + m+1)! 

Using the orthogonality relations for the Racah polynomials, see [21 p. 344], [T4"j §1.2], 

^2(m + n + l-2t) R k (X(t); 0, 0, -m -1,-n - 1) Ri{X{t); 0, 0, -m -1,-n- 1) 



t=o 



(n + l)(m + l)(m + 2) fc (n + 2) fc 
*'' (2fc + 1) (-m) fe (-n) fe 
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we find the following explicit expression for (3 k (m,n) 

1 (2k + 1) (-m) k (-n) k 



f3 k (m,n) 



C k {m, n) {n + l)(m + 1) (m + 2) k {n + 2) k 

(A.3) 

y^(m + n + 1 - 2t)R k (X(t); 0, 0, -m - 1, -n - l)a t (m, n)- 



m 
X 

t=0 



Now Proposition 12. 2[ and hence Theorem 12.11 follows from the following summation and 
simplifying the result. 

Lemma A.l. For £ e |N, n,m,k GN wi/i 0<k<m<nwe have 

X^" 1 ) r - ,J " tj -(m + n + 1 - 2t)R k (X(t); 0, 0, -m - 1, -n - 1) 



t=o 



(n + 2) m _ t t! 



(2^ + A; + l)!(2^-A;)! (n + 1) 



v ' (2£+l)! m!(2£-m)! 

Proof. Start with the left hand side and insert the 4.F3-series for the Racah polynomial and 
interchange summations to find 

y> (-A:) 3 - (fc + 1) 3 - (n-2i) m y> t (-l-Tl-m)t 

Z. j!j!( _ m ). ( _ n ). ( n + 2) m ^- l J (2^-n-m + l) t 



Tj ( _1 ) ("' 1 " " ' -' l\-i /.,»• 



-l)'(m + n + 1 - 2t)(-t)j(t -m-n-V 



Relabeling the inner sum t = j + p shows that the inner sum equals 

no (—l — n — m)oi , „ . . 
^ J (L-n-m + l) J ( " 2 '" lMl + m + W " 2i) 

'""' :-l-n-m + j) p (-2£-l+j) P 



x 



£ 



o (2£-n-m + l+j) p p! 



(1 + |(-1 - m - n + 2j)) P (-1 - m - n + 2j) p 



and the sum over p is a hypergeometric sum. Multiplying by (-_^w J _^ p the sum can be 



(\(-l-m-n + 2j)) p (-1 - m - n + j) p 

{j-m)p(j-n) p 
iij ouiii. ivj.uiLijJij'iiig, uy 

written as a very-well-poised 5 F 4 -series 

1 + | (— 1 — m — n + 2j) , — 1 — m — n + 2j, — 1 — 21 + j, j — m, j — n 

2 1 



|(— 1 — ?ri — n + 2j), 2£ — n — m + j + 1, j — n, j — m 



_ (-m-n + 2j) m _j (-m + 1 + 2£) m _j 
(-m - n + j + 1 + 2^) m _j (-m + j) m -j 

by the terminating Rogers- Dougall summation formula [U §4.4]. 
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Simplifying shows that the left hand side of the lemma is equal to the single sum 



fa-2*) TO + ^ {-n-m) m (2£ - I - m) 



(n + 2) m (2£ — n — m + l) m m\ 



k 



£ 



-fc)i (A; + 1),- {-2£ - 1) 



x - ■■ ^ ;j 



which can be summed by the Pfaff-Saalschiitz summation [2[ Thm. 2.2.6], [HI (1.4.5)]. This 
proves the lemma after some simplifications. □ 

Appendix B. Moments 

In this appendix we give an explicit sum for the generalised moments for W. By the explicit 
expression 

U r (x) = (r + 1) 2 F 1 I 3 ; — 
we find 



(1 - x) n U r (x)Vl - x 2 dx 

fr , n v brMr + ^ o-fcon+^ r^ + fc + 1) r (l) 

fr | n2 n +2 r(n+|)r(|) /-r,r + 2,n+§ 
(r + 1)2 r(n + 3) 3F2 v §,n + 3 '* 

= (r + 1)2 n +2 r (™ + D r (!) (-n)r 



r(n + 3) (n + 3) r 

using the beta-integral in the first equality and the Pfaff-Saalschiitz summation [2J Thm. 2.2.6], 
[Tlj (1.4.5)] in the last equality. For m < n, the explicit expression f 1 1.11) gives the following 
generalised moments 

- 1 ,T(p + §) r(|) (21 + 1) {2£ -m)\m\ 



J (l-x) p W(x) nm dx = 2 p+2 - 



T(p + 3) n + l (2f)! 

x £(-lH"-^(* + *-'>'(n + m - 2* + 1)±%^ (B1) 

f^ V ; (n + 2) m _ 4 t! v ; (p + 3) n+m _ 2< 
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